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Abstract. We prove that the sonic surface of axi- 
symmetric meridional stationary flows is always attached 
to the accretor, however small, if the adiabatic index of 
the gas is 7 = 5/3. 

Using local expansions near a point-like accretor, we ex- 
tend Bondi's classification of spherically symmetric flows 
to axisymmetric flows, introducing the possibility of angu- 
lar sectors reached by no flow lines, and singular directions 
of infinite mass flux, in addition to the angular regions of 
subsonic and supersonic accretion. For 7 < 5/3, we show 
the impossibility of subsonic accretion onto a point-like 
accretor when the entropy of the flow is not uniform. The 
special case 7 = 5/3 is treated separately. 
We analyse the influence of the adiabatic index and Mach 
number of the flow at infinity on the mass accretion rate 
of shocked spherical flows. We propose an interpolation 
formula for the mass accretion rate of axisymmetric flows 
as a function of the Mach number and the adiabatic index, 
in the range 9/7 < 7 < 5/3. 

Key words: Accretion, accretion disks - Hydrodynamics 
- Instabilities - Shock waves - Binaries: close - X-rays: 
stars 



1. Introduction 

Numerical simulations of the Bondi-Hoyle-Lyttleton 
(hereafter BHL) accretion flow, in 2-D and 3-D, have not 
only enabled the determination of the mean rate of accre- 
tion of mass, linear and angular momentum, but have also 
significantly modified our understanding of its dynamics. 
They have revealed that the accretion column foreseen by 
Hoyle & Lyttleton (1939) widens to form a detached bow 
shock if the adiabatic index of the accreted gas is 7 = 5/3 
(Hunt 1971) or 7 = 4/3 (Hunt 1979). The shock structure 
observed in simulations was used by Eadie et al. (1975) to 



interpret the X-ray observations of Vel X-1. 
Real situations like wind accretion in X-ray binaries re- 
quire taking into account the density and velocity gradi- 
ents in the incoming flow (see Ishii et al. 1993 and ref- 
erences therein), the effects of rotation (see Theuns et al. 
1996 and references therein) as well as the intricacies of 
radiative transfer (see an overview in Blondin 1994). For 
the sake of a better understanding, we restrict our present 
investigation to the axisymmetric meridional accretion of 
a gas with uniform adiabatic index onto a gravitating ac- 
cretor with constant velocity. 

Although this configuration is highly simplified, a com- 
pletely unexpected non-axisymmetric instability was dis- 
covered in numerical simulations for supersonic flows 
(Fryxell & Taam 1988). This instability could explain the 
observed variability of the radio or X-ray luminosity in 
some astrophysical systems (e.g. Ruffert & Melia 1994). 
However, its is not easy to establish the consistency of the 
results of the few published 3-D numerical simulations 
which depend on the details of the numerical method: 
no instability appeared in the 3-D SPH simulations by 
Boffin (1991) whereas it was observed in the 3-D Eule- 
rian simulations by Matsuda et al. (1992 and references 
therein) and Ruffert (1995 and references therein). Dome- 
shaped structures frequently appeared with Roe's method 
(approximate Riemann solver) of Matsuda et al. (1991), 
but not in the nested grid Piecewise Parabolic Method 
(PPM) of Ruffert (1991). Since numerical viscosity ap- 
pears to play an important role, a sufficient resolution of 
the simulation is crucial. In this respect, numerical simu- 
lations encounter the following difficulty: the smallness of 
the accretor in simulations is limited by the computational 
power, because the shortest timestep of the simulation is 
usually determined by the region of high velocities in the 
vicinity of the accretor. The ratio of scales between the ac- 
cretion radius and the radius r^, of a compact star, moving 
with a supersonic velocity v is typically of order: 
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where rschw. is the Schwarzschild radius. Note that in the 
case of accretion onto a neutron star, the relevant accretor 
boundary is the boundary of its magnetosphere, typicaUy 
a factor 10 — 100 larger than the Schwarzschild radius. 
Some of the most recent 3-D simulations (Ruffert 1992, 
1994a) could allow a ratio up to rA/r* ~ 10^ (the full 
domain of the simulation extends up to ~ 16 accretion 
radii), by using a PPM code with multiple nested grids. 
Such an accreting sphere is still more than 100 times larger 
than rgchw., and typically 10 times larger than the mag- 
netosphere of a neutron star. One would like to be able to 
predict, on an analytical basis, which changes are likely 
to occur if the size of the accretor is further decreased. 
The role of the boundary conditions at the surface of the 
accretor is closely related to the position of the sonic sur- 
face, and we shall use analytical arguments to obtain some 
insights on this question. 

An analytical description of axisymmetric stationary 
flows is also a first step towards the analysis of their sta- 
bility, which we postpone to a future paper. Bisnovatyi- 
Kogan et al. (1979) found self-similar stationary solutions, 
and one would like to know to what extent any generic 
stationary solution would behave in the same way as their 
solution. More generality is gained by studying, in the 
same spirit as Theuns & David (1992) for the spherical 
case, the first order behaviour of all quantities near the 
accretor for an axisymmetric flow, and in particular their 
departure from sphericity. 

Some analytic formulae have been proposed to fit the 
mass accretion rate observed in simulations (Hunt 1979; 
Ruffert 1994b), but none of them filled the gap between 
the well studied cases of axisymmetric accretion of dust 
(Hoyle & Lyttleton 1939; Bondi & Hoyle 1944) on the one 
hand, and spherical accretion of a gas (Bondi 1952) on the 
other, for an arbitrary value of the adiabatic index . Con- 
tinuing in this direction, we add to these reference models 
the spherical accretion of a shocked gas with uniform adi- 
abatic index, in order to see the effect of both the kinetic 
and the thermal energies, and extend this approach to ax- 
isymmetric meridional stationary flows. 

We first recall in Sect. || the general equations gover- 
ning a stationary fiow and formulate the axisymmetric 
problem with a single equation in cylindrical coordinates. 
The shape of the sonic surface is analysed in Sect. |^. We 
establish some properties of the fiow by performing local 
expansions in the vicinity of the accretor in Sect. ]^ We 
present in Sect. || the analytical solution of the spherical 
accretion from a gas with constant kinetic and thermal en- 
ergies at infinity. An interpolation formula for the mass ac- 
cretion rate of axisymmetric meridional fiows is proposed 
in Sect. ^. We summarize our conclusions in Sect. |^. 

2. Equations of a stationary flow in a gravitational 
potential 



2.1. General equations 

The continuity equation for a gas of velocity v and density 
p is: 



V • (pv) = . 



(2) 



We restrict our study to the case of a gas of uniform adi- 
abatic index 1 < 7 < 5/3. The entropy S* of a gas with 
pressure P is defined as: 
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(7 - 1)/^ 



log- 
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in which TZ is the gas constant and /i is the mean molecular 
weight. The entropy is constant along each flow line {v ■ 
VS" = 0), except when matter passes through a shock. 
The temperature T and the square of the sound velocity 
are simply related through a multiplicative constant: 



7P _ -fTZT 
P 



(4) 



For the sake of simplicity, we shall set the ratio TZ/p, = 
1 everywhere in what follows. Since the mean molecular 
weight is constant throughout our study, this convention 
is equivalent to changing the units of temperature and 
entropy, without loosing any generality. 
Using Eq. (|3|), the pressure force can be written as: 
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The Euler equation for a stationary flow with velocity v 
around an accretor of mass M is then: 



V— (V X t>) X t> = V 

2 r 7 — 1 



Vc^+TVS, (6) 



in which r is the distance to the centre of the accretor and 
G is the gravitational constant. 

The Bernoulli equation is obtained by integrating the Eu- 
ler equation along a fiow line. It states that the quantity 
B 



B 



GM 



7-1 



(7) 



is constant along each fiow line. The three terms of this 
sum are respectively the kinetic, the thermal, and the 
gravitational energies. 

The vorticity of the fiow is simply defined as 



to = V X t> 



(8) 



By using the definition of the Bernoulli constant (Q) in the 
Euler equation (^, we obtain the following equation for 
any stationary flow, 



WKv = TVS - VB . 



(9) 
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This formula is particularly useful because it expresses the 
vorticity in terms of invariant quantities B and S. 
The Bernoulli constant B remains constant along a flow 
line, even through a shock. By contrast, the entropy S in- 
creases through a shock. 

In the classic BHL case of a flow with constant and uni- 
form velocity at infinity, the Bernoulli constant is uni- 
formly constant (VB = o), and so does not contribute to 
the vorticity: 



wxv = TVS 



(10) 



If, in addition, the flow is isentropic at infinity, the only 
possible source of vorticity gradients is the shock front. 
Note that Eq. (|l^) implies the conservation of entropy {v ■ 
VS* = 0). Thus the only remaining differential equations 
are the continuity Eq. (||) and the vorticity Eq. (|l0|). 

2.2. Third-order partial differential system of an axi- 
symmetric stationary flow in cylindrical coordinates 

We consider the axisymmetric stationary flow of a gas 
with uniform density and constant velocity v^o \\ x with 
respect to the accretor. The geometry of the flow lines is 
described by a single function, the angle f3{r,9) between 
the radial direction and the flow line, defined by: 



tan/3 : 



We 

Vr 



(11) 



We define by n the direction perpendicular to the velocity, 
in the plane containing the symmetry axis. Denoting by 
er,ee the unit vectors in cylindrical coordinates, we can 
write: 



V = —v{cos Per + sin f3ee) , 
n = sin f3er — cos Peg . 



(12) 
(13) 



When useful throughout the paper, we denote by f the di- 
mensionless distance normalized to the modified accretion 
radius te, defined as: 



GM 



GM 



Boo ^ 



(14) 



It coincides with the classic accretion radius ta for large 
Mach numbers, when the kinetic energy largely exceeds 
the thermal energy: 



2GM 



(15) 



We shall denote distances normalized to te by a bar. 
Using Eqs. (||) and (^, all relevant fiow quantities can be 
rewritten as functions of S", A^, /?: 



V = -M 



2Boo(l + r) 
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(cos/3er + sin /See) , (16) 



2goo(l + f) 

2Boo(l + r) 



P = e 



2^0^(1 + f) 



The azimuthal component of the Euler Eq. 
by r, is: 



drv0 
dr 



dvff 



1 dP 



de pde 
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(17) 
(18) 
(19) 

), multiplied 
(20) 



We denote by w the third and only non-zero component 
of the vorticity vector: 



w 



1 drvg 
r dr 



1 dvr 
r~de 



(21) 



By combining the vectorial vorticity Eq. ( |l0| ) and the 
scalar continuity Eq. (H), we obtain a third order differ- 
ential system in terms of /3, A^^, S only, with variables r, 
(see Appendix ^) . This system is hyperbolic in regions 
where A4 > 1, and elliptic where < 1. It requires three 
functions P^{6), Ml{9), S^,{9) as boundary conditions on 
the surface of the accretor. It should be noted that if the 
entropy is uniform, the fiow is potential {v = V$v), and 
the Bernoulli and the continuity equations provide a sec- 
ond order partial differential equation on <&v(?', 9). In this 
case, the boundary conditions at the surface of the accre- 
tor consist of two functions of 6, namely $v(?' = ^) and 
its radial derivative. 

2.3. Further reduction of the differential system 

We can simplify the system by using the modified stream 
function introduced by Crocco (1936), defined as follows: 



a* 

'09 
dlogr 



= (2GM) 2(T-i)c— r^'y^sine < 
9^- 



tan/3 



d9 



(22) 
(23) 



Where there is no shock, 'J is conserved along fiow lines 
(because v ■ V^* = 0). It is also conserved across shocks, 
since the derivatives of ^' have only a finite discontinu- 
ity there. The entropy can be expressed as a function 
of 4*, and thus becomes a function of a single variable 
S{r,9) — S{'^{r,9)), defined once and for all along the 
shock. 

Let us index the flow lines with their distance w to the 
symmetry axis (i.e. their impact parameter) at infinity. 
Only the gas within a cylinder of radius zuq is accreted, 
so that the total mass accretion rate is M — irvo'^pooVoo. 
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At infinity, the mass flux inside the cyhnder of radius m is 
directly proportional to ^'(nj) through the following prod- 
uct: 



(2GM)- 

7 



^'(n7)e 



(24) 



For the sake of clarity of the equations, let us define: 
9* 



e^, — r 



dr 



(25) 



Because the total mass accretion rate is finite, we deduce 
from Eq. (^4|) that ^ is bounded in the vicinity of the 
accretor, and linv^o e* = 0. 

Using the definition of '5 and the Bernoulli equation, we 
can formally write the Mach number M. as a, function of 
the gradient of through the following algebraic equation: 



one of the two roots into Eq. (pq), we obtain an implicit 
differential equation of second order containing explicitely 
only ^1 and its derivatives of first and second order. 
This compact formulation will allow us to perform in 
Sect. ^ a direct analysis of the behaviour of the solution 
in the vicinity of r = 0. 

3. The sonic surface of axisymmetric stationary 
flows 

3.1. Three topologies of axisymmetric accretion 
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At a given position (r, 9) , the maximum of the right hand 
side term is reached at = 1. This defines a maximum 
value for the gradient of the stream function, i.e. for the 
local mass flux according to Eq. ( |2^ ) . For any value of V^* 
below this maximum, two Mach numbers are solutions of 
the Eq. (p6|). As for spherically symmetric flows classified 
by Bondi (1952), one solution is subsonic, while the other 
is supersonic. 



Equation (^6|) can be derived with respect to r or to 
express the partial derivatives of M with respect to r or 
as function of A4 and the derivatives of The vorticity 
Eq. (^0|) can therefore be expressed with derivatives of 5* 
only, as follows: 



+ f 



df 
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2"^ 2(7-l)(l-A^2) 

2 + (7-l)A^2 



1 + f 



2(7-l)(l 
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de de 



9* 



86 J 



C(27) 



This equation can be viewed as an explicit partial differ- 
ential equation of second order on ^E", with Ai implicitly 
depending on V^* through the algebraic Eq. (p6|). In an 
equivalent way, Eq. (27) can also be considered as a second 
order polynom in A^2^ whose two roots explicitly deter- 
mine 7V(2 as a function of 'J and its derivatives. Inserting 



Fig. 1. Three topologies of axisymmetric accretion. The sonic 
surface is represented as a dotted line. It touches the accretor at 
an angle Oso with the symmetry axis for the types SF and FSF. 
For the type FSF, the solid line represents the shock surface 
reaching the accretor at an angle 9sh with the symmetry axis. 



Numerical simulations of axisymmetric accretion with dif- 
ferent Mach numbers, adiabatic indices and accretor sizes 
have exhibited three different topologies of the sonic sur- 
face. These three topologies are illustrated in Fig. |l|. 

- Type F (fast accretion) : the sonic surface is detached 
from the accretor. This occurs for subsonic fiows for small 
enough accretors (Ruffert 1994b, 1995, 1996). 

- Type SF (mixed slow/fast accretion): the sonic sur- 
face reaches the surface of the accretor. It subtends an 
angle Ogo relative to the axis of symmetry at the position 
of the accretor. This is the case for subsonic flows with 
large accretors, and supersonic flows with 7 = 5/3 (Mat- 
suda et al. 1992; Ruffert & Arnett 1994; Ruffert 1994a). 

- Type FSF (mixed fast /slow/fast accretion): both the 
sonic surface and the shock surface touch the surface of 
the accretor. They respectively subtend the angles Sgo and 
9sh relative to the axis of symmetry. This occurs for all su- 
personic flows when the accretor is large enough, for all 
the known supersonic 3-D simulations with 7 close to one 
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(Ishii ct al. 1993; Ruffert 1996). More complicated topolo- 
gies were observed by Hunt (1979) and Ruffert (1995) in 
simulations of supersonic flows with 7 = 4/3, where the 
bow shock is linked to the accretor by a conical shock. 
These can also be classified as Type FSF, with the differ- 
ence that the entropy in the supersonic flow ahead of the 
attached shock is not uniform. 

The topology of the flow in numerical simulations appears 
to depend strongly on the size of the accretor: all the 
supersonic flows are of Type FSF if the accretor is big 
enough, but some become Type SF or F when the accre- 
tor size is decreased. In order to spare computing time, 
which increases with decreasing accretor size, one would 
like to use the largest accretor size leading to a physically 
relevant accretion process (cf. Eq.||). Type F flows are nu- 
merically easier to handle, because unphysical effects due 
to a possible reflection of waves at the boundary of the 
too large accretor can be avoided. 

3.2. The shape of the flow tubes at the sonic radius 

Let us consider a flow tube of cross section $, 
parametrized along the tube by the curvilinear variable 
We denote by AS the entropy that this flow tube gains 
passing through a shock, and by poo and Coo the density 
and the sound speed at large distances from the accretor, 
before the shock. The density p and sound speed c after 
the shock satisfy the following equation: 



-AS 



/Coo\ ^ _P_ 
\ C J Poo 



(28) 



Using the Bernoulli equation (Q) and the entropy equation 
(p8|), we can write the momentum as a function of the 
Mach number M. and the position r: 



pv 



Mpoo AS 



(7 - l)Boo (1 



1 



(29) 



The continuity equation (^ states that the mass flux 
is conserved along the flow tube: 



d_ 
dl 



{pv<^) = . 



By substituting Eq. 
given flow line: 



91og$ 
dlogr 



7 + 1 



(30) 

]) into Eq. (pO[), we obtain, along a 
1-X2 9iogX2 



2(7-l)l-Hf 2 + (7-l)A^2 aiogr 



,(31) 



where the derivatives are calculated along a flow line, i.e. 
d/dr = {d/dl){dl/dr). The regularity of dM/dl at the 
sonic point A^(fso) = 1 implies: 



51og$ 
dlogr 



= 2-f 2 



(32) 



Fig. 2. The position of the sonic point depends on the shape of 
the flow tube: for an axisymmetric flow tube, it is shifted from 
the spherical position depending on whether its cross section 
increases faster or slower than the square of the distance to the 
accretor 



We have denoted by rg the radius of the sonic sphere in 
the spherical case: 



37 



GM 



(33) 



It is remarkable that the spherical sonic radius is a func- 
tion of the total energy B^o only, and is independent of 
the entropy of the flow. 

The regularity condition (^2|) is a constraint on the shape 
of the flux tube at the sonic radius. Note that this con- 
dition is valid for any stationary flow, with no hypothesis 
on its symmetry. Let us illustrate it for an inflowing ax- 
isymmetric flow tube: the spherical case corresponds to 
radial flow lines, and thus to a conical flow tube. If the 
sonic radius is smaller than the spherical one, the cross 
section of the flow tube increases faster than the square 
of the distance to the accretor, thus the azimuthal com- 
ponent vg{r so) of the velocity at the sonic radius must 
be oriented towards the axis of symmetry. Conversely, if 
the sonic radius is larger than the spherical one, the cross 
section of the flow tube increases more slowly than the 
square of the distance to the accretor, and we(fso) must 
be oriented away from the axis of symmetry (see Fig. H). 

3.3. A geometrical property of Type F flows 

We prove in Appendix ^ the following property V: 
If the sonic surface is detached from the accretor, it must 
intersect at least once the sonic sphere of radius tq of the 
corresponding spherically symmetry flow with same en- 
ergy. 
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face and the shock surface appear to be attached to the 
accretor. Although necessary, the condition r^, < tq is far 
from being sufficient. 



3.3.1. Consequence for flows with 7 ~ 5/3 

Since ro(7 = 5/3) = 0, the sonic surface of flows with 
7 = 5/3 must always reach the surface of the accretor. 
Following the nomenclature introduced in Sect. 3.1, the 
topology of the accretion with 7 = 5/3 is either of Type 
SF or FSF, for any accretor size and any Mach number 
at infinity. 



Fig. 3. If the sonic surface (dashed line) of the axisymmetric 
flow is detached from the accretor, it must intersect the sonic 
sphere (dotted line) of the equivalent spherically symmetric 
flow. 



ro is defined in Sect. 3.2 by Eq, 



33). This property can be 
understood intuitively from Fig.T?, by using the geomet- 
rical interpretation of Eq. (p2|) in Sect. 3.2, in the simple 



case where the azimuthal velocity is of constant sign along 
the sonic surface in the interval < 6 < n.li the azimuthal 
velocity is directed away from the axis on one side of the 
accretor (for example, 7'so(0) > ro), it will be oriented 
towards the axis on the other side (rso(7r) < ro), thus im- 
plying by continuity an intersection of the sonic surface 
with the sonic sphere. We show in Appendix ^ that this 
property is true for any axisymmetric flow, without any 
hypothesis on the distribution of azimuthal velocities. 
Former studies of the BHL accretion have dealt with two 
set of length scales, one associated to the Hoyle-Lyttleton 
accretion of dust (the accretion radius ta), and the other 
associated to the spherical accretion (the Bondi radius re 
and the corresponding sonic radius ro(7Woo = 0))- None of 
these scales took into account both the kinetic energy and 
the thermal energy of the gas. In particular, the signifi- 
cance of the sonic radius ro(A4oo = 0) for axisymmetric 
flows has always been dubious. The property V establishes 
the fundamental meaning of the spherical radius 7-0(^^00), 
in relation with the sonic surface, for axisymmetric flows 
of Type F. 

The geometrical property V is useful for numerical simula- 
tions of the BHL flow, in order to fix the maximum size of 



the totally absorbing accretor. As discussed in Sect. 3.1 
since the sonic surface must cross the sonic sphere, the 
accretor size should at least be smaller than the spher- 
ical sonic radius tq. However, although the sonic sphere is 
larger than the accretor size used in numerical simulations 
with 7 = 4/3 and 7 = 1.01 (see Ruffert...), the sonic sur- 



3.3.2. Consequence for nearly isothermal supersonic fiows 

The Hoyle-Lyttleton limit for the mass accretion rate at 
high Mach numbers requires that the accreted matter 
comes from a cylinder with a radius scaling like the ac- 
cretion radius (Hoyle & Lyttleton 1939). If 7 is close to 
unity, the sonic spherical radius ro is A^^/4 larger than 
the accretion radius rA: 



ro_ 

rA 



5-37 



M 



Mi 



for 7 = 1 



(34) 



If such a flow were of Type F, the geometrical property 
V would imply that the sonic surface goes out as far as 
A4'^/A times the accretion radius. This leaves only two 
possibilities for a stationary flow: 

(i) the accretor is able to keep material bound up to 
distances comparable to ro ^ rA, and the mass accre- 
tion rate is inflnitely larger than the Hoyle-Lyttleton limit. 
This is the case of the spherical flows we shall study in 
Sect. |. 

(ii) the accretor cannot retain material at such large 
distances, and the flow is of Type SF or FSF, as observed 
in numerical simulations with 7 close to one by Ishii et al. 
(1993) and Ruffert (1996). 

Let us note that the artificially constructed Type F geom- 
etry, where the distance of the stagnation point exceeds 
ro, while the accreted matter collects within a cylinder of 
radius ~ rA, would be strongly unstable to the Kelvin- 
Helmholtz instability: the growth rate would scale like 
Coo/rA while the time needed to reach the accretor 
would scale like ~ Coo/rQ. 

Consequently, we rule out the possibility that the sonic 
surface is detached in the case of an isothermal stationary 
flow. According to Eq. (34), the spherical sonic radius is 
larger than the accretion radius at high Mach numbers if 
(see also Fig. ||) 



1 < 7 < 



9 



(35) 
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This angle can be smaller than tt. In particular, 9c < 
9sh for flows of Type FSF. Without loosing any signif- 
icant physical generality, we can assume that ]O,0c[C 
{6io(]0,roo[)} in the BHL flow. 



Fig. 5. Illustration of our notations 6c,vj ,6o['co) 



Fig. 4. Comparison of the accretion radius ta with the sonic 
radius ro depending on the values of 7 and Moc- The thin 
solid line indicates the region of parameters where the kinetic 
energy is equal to the thermal energy. The value of rA/ro is 
indicated on the dashed and dotted curves 

4. Asymptotic expansions near the accretor 

Spherically symmetric accretion flows with 1 < 7 < 5/3 
are separated into two distinct classes, depending on 
whether the flow is subsonic or supersonic at the accre- 
tor boundary (Bondi 1952). We learnt from Eq. ( p6|) in 
Sect. 2.3 that these two branches of solutions also exist for 



axisymmetric flows, for a given value of the local mass flux 
at the surface of the accretor. The azimuthal parameter 0, 
however, introduces an additional degree of freedom, such 
that the solution might be subsonic in an angular sector, 
and supersonic in an other (i.e. for a sonic surface attached 
to the accretor). We must also distinguish between angu- 
lar directions towards which the flow lines converge (direc- 
tion of "accreting" flow lines), and directions along which 
no gas is accreted (direction of "passing" flow lines). By 
conducting a local analysis, we aim to clarify the various 
possible configurations. 

^.1. Definition of regular and singular accretion 

For each flow line (indexed by w) reaching the accretor, we 
denote by 9q{tu) the angle between the velocity vector v 
and the axis of symmetry, when r 0. Our notations are 
illustrated in Fig. ^. We define the direction 6* = along 
the axis of symmetry on the downstream side of the accre- 
tor, so that the function 9q{izi) decreases monotonically in 
the interval [0,i37o], with 6*0(0) — tt and 610(^370) = 0. We 
define 9c as the largest azimuthal angle, at r — s- 0, among 
all the converging fiow lines, apart from the axis of sym- 
metry: 



9c=M&x{9o{]Q,mo[)} 



(36) 



Near the surface of the accretor, let us denote by 
27rTOo(6') sin 6'd0 the local fiux of mass between the az- 
imuthal angles 9 and 9 + d9: 



moi9) 



lim - 



lim - 



-r pVr , 



(2GM)- 



1 



sin0 89 



lim 



e cos (3 
5-37 



7^-1 7'2(-y — 1) 



2GM{\ + f) 



2(7-1) 



2 

7-1 



(37) 
(38) 

(39) 



Since the stream function is bounded in the flow (Eq. 
integrating Eq. (23) with respect to r implies: 

,9* 



lim /3-„ 

r,-*0^ d9 



for 6* e [0, tt] 



(40) 



Thus, if the mass accretion rate mo is locally not zero, the 
flow lines become locally radial when r — s- (no "spiral- 
ing" of the flow lines). On the other hand, if the flow lines 
are not radial, Eq. (38) implies that the mass accretion 
rate is zero in this direction (passing flow lines). 
Note that along a direction of accretion, we can consider 
a flux tube of cross section around the flow line such 
that 9q{w) — 9, and write Eq. (37) as follows: 

^2 



TOo(6'o(n7)) 



'(pt;$ro) lim -— 



(41) 



The quantity niQ is consequently finite only if the fiow 
tubes are asymptotically conical {^^{r) ^ r^). This leads 
us to call regular the accretion onto a point-like accretor 
such that Too remains finite, and singular when it is locally 
infinite. 

We can also express the local mass accretion rate in terms 
of the mass flux between the cylinders of radii tu and w + 
dm at infinity: 



TOo(6'o(n7)) 



sin 9o dOp 
d-co 



(42) 
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Thus singular accretion corresponds to a localized az- 
imuthal angle 9o where d6'o/dn7 = 0. An example of 
singular accretion is the self-similar solution found by 
Bisnovatyi-Kogan et al. (1979), where the whole accre- 
tion is along a single azimuthal angle 0k = 0o(]O, wo])- 
For a regular flow (^^(r) ^ r^), we take the limit of 
Eq. (^ij) when r ^ and obtain, along a given flow line. 



dlogM' 



5 - 37 2 + (7-l)7W2 



Slogr 2(7-1) 1-A^2 



(43) 



This equation shows that for 7 < 5/3, the Mach number 
along regular flow lines either tends to zero or to infinity 
when approaching the accretor. As in the case of spherical 
accretion (Bondi 1952) no intermediate finite Mach num- 
ber is possible. Moreover, this equation proves, by conti- 
nuity, that a mixed subsonic/supersonic accretion requires 
a singular accretion at the transition if 7 < 5/3. The sonic 
surface may still reach the accretor in a regular flow, by 
separating a region of supersonic accretion from a subsonic 
region of passing flow lines, as we shall see in Sect. 4.5. 
We will now examine the different regions of passing flow 
lines, supersonic and subsonic regular accretion succes- 
sively, for 7 < 5/3, and treat the case 7 = 5/3 separately. 
We shall use the generic name e(r) for a function which 
tends to zero with r, and the subscript for the limit of 
functions when r 0. 

4-2. Region of passing flow lines 9 £\9c,tt[ for 7 < 5/3 

Since the flux of mass is zero, Poi9) cannot vanish in 
the interval ]0c,7r[. According to our sign convention in 
Sect. 4.3.3 and the definition of f3 in Eq. (|ll|), 



f3o{9) > ior9e]9c,7T[ 



(44) 



Since Po{Tr) ~ P{9c) = 0, the function /3o must have a 
maximum in the interval ]9c, 7r[. 

(i) If the shock is detached, the entropy is finite in the 
flow. According to Eq. (39), the Mach number must either 
diverge or tend to zero in order to satisfy — 0. 
Let us suppose that Al — > 00. We use Eq. (16) and obtain: 



z« = -^sin/3ofl + 2^)(l + .(r)). 



2r 



89 



(45) 



Using the vorticity equation ( [lO| ) , the entropy gradient can 
be expressed as: 

V5 = sin/3o (1 + 2^) (1 + e(r)) . (46) 

Integrating this equation over a segment between 9 = tt 
and 9 = ^max, for r — > 0, contradicts the fact that the 
entropy is bounded for a flow with a detached shock. Con- 
sequently the Mach number tends to when r — > 0, and 



^so < 9c. We conclude that all the flows with a detached 
shock and 9c < tt are of Type SF. 

(ii) If the shock is attached, the incoming supersonic 
velocity diverges near the accretor like r^^/^ when r — > 
0. We decompose it into a parallel and a perpendicular 
component with respect to the shock. The Mach number 
M2i_ associated to the component of the velocity, after 
and perpendicular to the shock, is subsonic according to 
the Rankine-Hugoniot jump conditions (see e.g. Landau 
& Lifschitz 1987): 



lim M2 



7shj = 



7-1 
27 



(47) 



The Mach number M2\\ associated to the component par- 
allel to the shock can be supersonic, but not arbitrarily 
large: 



limTWalK^'sh) 



7 + 1 



1 



[27(7 - 1)]2 tan/3(6'sh) 



(48) 



The total Mach number immediately afler the attached 
shock is consequently finite when r 0. This result will 
be important in order to study the local stability of such 
flows to the Rayleigh-Taylor instability. 

4-3. Region of regular supersonic accretion for 7 < 5/3 

If — > 00 and toq is finite, Eqs. (16)-(19) can be ex- 
panded as follows: 



f2GM\ ^ , 
Vr = - \ j (l + e(r)) 



p = TUq 



/ 1 



72 



(2GMr3 
1 



— (1 + ^0) 



p = e(--i)^«m2 ( ) (l + e(r)) 



M 



1-1 

So\ — 



mo 

00 when r — > 



2GMr3 
(2GM)^ , , . , 

'Y 2 7" 4 



(49) 
(50) 
(51) 

(52) 
(53) 



The first order azimuthal velocity can be estimated using 
Eq. (|o|) in the azimuthal Euler Eq. (pO|). If the pressure 
is not spherically symmetric to first order. 



1 dP 



pd9 2 

/ 4-3t\ 

ve ^ O (r— j 



(54) 
(55) 



Vr being negative, the azimuthal velocity is oriented to- 
wards the region of increasing pressure. The azimuthal 
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pressure force is balanced by this dynamical term. We de- 
duce from Eqs. ( pi] ) and (55) an upper bound for the de- 
parture from spherical symmetry of the function (3: 



lim 

r— tO 



91og/3 
i91ogr 



> 



37 



(56) 



By studying the leading terms of Eq. (27) when r ^ 0, 
we derive in Appendix |c| a lower bound for the departure 
from spherical symmetry of the function /?, if the entropy 
of the flow is not uniform: 



lim 



d\ogr 



<5-37 ifS"(*)7^0 



(57) 



On the other hand, the function [3 may approach much 
more rapidly than ^ r^-^^ if the entropy is uniform. To 
phrase it more qualitatively, the presence of an entropy 
gradient sets a lower bound on the "bending" of the flow 
lines near the accretor. If a parallel can be drawn with 
the simpler case of two-dimensional supersonic flow past a 
concave profile (e.g. Landau & Lifshitz 1987), this "bend- 
ing" could be related to the formation of long-lived ra- 
dial shocks in the supersonic region of flows with a non- 
uniform entropy, as observed in numerical simulations 
with 7 = 4/3 in 2-D (Hunt 1979) and 3-D (Ruffert 1995). 

Because the radial velocity (49) is spherically symmet- 
ric to first order, the azimuthal dependence of the accreted 
momentum pv comes mainly from the azimuthal depen- 
dence of the density p, which we decompose into the con- 
tributions of pressure and entropy, 



p(r, 9) — Pi ir, 0) exp ( — 



7-1 



S{r,B) 



(58) 



Thus we expect that for a flow with uniform entropy (no 
shock), the mass is mainly accreted from the hemisphere 
where the pressure is the highest. The geometrical inter- 
pretation of Eq. (|3|) implies that near the axis of sym- 
metry (i.e. for ~ and ~ tt), the azimuthal speed is 
oriented away from the axis where Vgo > tq, and towards 
the axis where Tso < tq-U one assumes that the sign of the 
azimuthal velocity does not change in the supersonic re- 
gion between the sonic surface and the accretor, Eq. (54) 
indicates that the highest pressure is on the side of the 
accretor where rgo > tq. 

For an isentropic subsonic flow, the natural orientation 
of the azimuthal velocity is set by the sign of the angu- 
lar momentum of the fluid element at infinity. A subsonic 
medium would flow from the left of the accretor in Fig. ^ 
Therefore the sonic surface is displaced opposite to the 
direction of the flow at infinity, and the maximum flux 
of mass is expected on the back hemisphere of the ac- 
cretor. Indeed, in the subsonic numerical simulations of a 
flow with 7 = 4/3 by Hunt (1979) in 2-D, and by Ruffert 
(1995) in 3-D, the sign of the accreted linear momentum 
is negative for the smallest accretor considered. 



If the flow is supersonic at infinity, the position of the 
shock along the axis is displaced in the direction of the 
flow, compared to its position for a spherical flow with 
same kinetic and thermal energies at infinity, because 
there is now a leakage of gas away from the accretor and 
the pressure ahead of the accretor is lower than in the 
corresponding spherical case. The sonic point along the 
axis is displaced accordingly. The sonic surface is conse- 
quently shifted in the direction of the flow, leading to a 
highest pressure on the front hemisphere of the accretor. 
The simplest configuration is the one in Fig. ^, where the 
supersonic flow now comes from the right of the accretor. 
Nevertheless, the entropy generated by the shock de- 
creases with increasing distance from the axis of symmetry 
(see Appendix ^ , so that entropy and pressure contribute 
in opposite directions to the density in Eq. (p8|), preclud- 
ing a direct conclusion about its azimuthal distribution 
near the accretor. However, the accreted linear momen- 
tum is negative in the numerical simulations by Ruffert 
(1996 and references therein) indicating that the entropy 
effect is stronger than the pressure effect. 

4- 4- Region of regular subsonic accretion for 7 < 5/3 

If ^ and toq is finite, Eqs. (16)-(19) can be written 
as follows: 



Vr 



-So 



7 



(7 - l)GM 
(7 - 1)GM 



jr 

(7 - l)GM' 



3-27 

r 7-1 (1 -I- e(r)) 



(1+6(0) 



(l + e(r)) 



P 



(7 - l)GM 



jr 



(l + 6(r)) 



M = 2(7-1) (l-He(r)) 

[{j-l)GM]~ 

— > when r . 



(59) 
(60) 
(61) 
(62) 

(63) 



We use Eqs. (59) and ( [40| ) to estimate the dynamical 
terms balancing the azimuthal pressure force in Eq. (20) 
and obtain: 



1 dP 



/ 2(3-2t) \ 



(64) 



Comparing this to the first order of P/ p ^ c^ given by 
Eq. (61), we conclude for 7 < 5/3 that the pressure in the 
regular subsonic region is spherically symmetric to first 
order when r —^ 0. 

Thus Eq. (62) implies that the subsonic flow for 7 < 5/3, 
with its first order expansions (59)-(63), is possible only 
if the entropy is uniform. 

In particular, the spherical solution proposed by Theuns & 
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David (1992), connecting a supersonic solution to a sub- 
sonic solution through a shock, cannot resist the slight- 
est asphericity if the accretor is point-like. The azimuthal 
pressure force due to the convergence of flow lines of dif- 
ferent entropies is too big to be balanced by dynamical 
forces in the regular subsonic region close to the accretor. 
Nevertheless, this constraint of the azimuthal Euler equa- 
tion is less stringent for large accretor sizes. One can also 
argue that the radiative processes near the surface of the 
star may homogenize the entropy, and allow for a spheri- 
cally symmetric subsonic accretion. 

The entropy in flows of Types SF and FSF is a strictly 
decreasing quantity along the shock: 

(i) In Type FSF flows, the velocity diverges to infinity 
when Tsh 0, thus the entropy created along the shock 
decreases strictly with distance from the accretor for Type 
FSF flows, precluding a subsonic accretion with a uniform 
entropy. 

(ii) For Type SF flows, we show in Appendix ^ that 
the entropy created along the shock strictly decreases with 
distance from the axis of symmetry. 

We conclude that there cannot be a region of regular sub- 
sonic accretion in shocked flows with 7 < 5/3, for a point- 
like accretor. The subsonic region may reach the surface 
of the accretor for any finite size, but the mass accreted in 
the shocked regular subsonic region 9 > 9so tends to zero 
when the accretor size decreases to zero. 

4-. 5. The sonic surface for Type SF and FSF shocked flows 
7 < 5/3 

Since no regular subsonic accretion is possible with a non- 
uniform entropy, all the regular flow lines are ultimately 
supersonic when they reach the point-like accretor {0c < 
9so)- Note that in the singular solutions by Bisnovatyi- 
Kogan et al. (1979), the Mach number tends to infinity 
along each flow line. 



the local mass accretion rate rriQ is always finite, its max- 
imum value being reached for A4o ^ 1: 



Together with the results of Sect. 4.2, we conclude that 
for a flow with a detached shock, and 7 < 5/3, 

0c = 0,0 ■ (65) 

According to Eq. (p^), the cross section of the flow tube 
reaching the sound speed at the surface of the accretor 
varies like: 



d log $ dl 



7-1-1 



dl d\ogr 2(7-1) ^ ^ ■ ^^^^ 

Comparing this equation to Eq. (pl|), we conclude that 
the azimuthal angle 0so cannot be a direction of regular 
accretion. However, 0so is not necessarily a direction of 
singular accretion. A regular flow simply requires that 0c ^ 
{0n(]O,n7o[)}. 

4-6. The particular case 7 = 5/3 

We denote by A^o(^) the limit of the Mach number along 
a converging flow line. For 9 g]0, 0c[, Eq. (39) implies that 



Too 



(Ml + 



< '5 



3\* G^M^ 



(67) 
(68) 



The pressure forces for 7 = 5/3 are strong enough to im- 
pede an arbitrarily high mass flux, even locally. Thus ac- 
cretion flows with 7 = 5/3 are always regular, and the 
maximum mass accretion rate is the spherical Bondi value. 
This result is particularly important to understand the in- 
stability observed in numerical simulations. Since the most 
unstable simulations were observed with 7 = 5/3, we can 
restrict our future stability analysis to the simpler case of 
regular flows. 

Using Eq. (16) and ( |4C| ) in the azimuthal Euler Eq. (pO[), 
we find that the azimuthal pressure force is bounded by 
Eq. (p4|). Comparing it to Eqs. (17) and (19), we conclude 
that the pressure is spherically symmetric to first order for 
£]0, 0c[- According to Eq. (p8[), the density is spherically 
symmetric to first order if there is no shock (uniform en- 
tropy) , but increases towards the back hemisphere if there 
is a shock. 

With Eqs. (16) to (19), we conclude that: 

(i) in the absence of a shock, the pressure, density, tem- 
perature, velocity and Mach number are spherically sym- 
metric to first order when r — s- 0. The accreted momentum 
is therefore spherically symmetric as well. This is con- 
firmed by the numerical simulations by Ruffert (1994b), 
with Aioo — 0.6 and 7 = 5/3. The accreted linear mo- 
mentum decreases to zero with the size of the accretor 
(see Fig. 2, model "SL", "SM", "SS"). Moreover, the uni- 
form Mach number A4o must be lower than or equal to 
one because of the property V. 

(ii) if a shock is present, the pressure is spherically 
symmetric to first order when r ^ 0, the temperature 
and entropy increase towards the upstream hemisphere, 
and the density, the velocity and the Mach number in- 
crease towards the downstream hemisphere. Consequently, 
more mass is accreted from the downstream hemisphere 
than from the upstream hemisphere. This is also confirmed 
by the numerical simulations with 7 = 5/3 by Ruffert 
(1994b), at AdoD = 1-4. This trend persists in the unstable 
cases, simulated by Ruffert & Arnett (1994) at A4oo = 3, 
and Ruffert (1994b), at Moo = 10. 

Thus we disagree with Hunt (1971), who chose the Mach 
number at the surface of the accretor tVJq in order to max- 
imize the mass accretion rate. As can be checked from 
Eq. (67), the highest mass accretion rate is obtained for 
A4o = 1. The Mach number at the surface of the accretor 
is, however, not a free parameter, but is constrained by 
the distribution of entropy. Let Si be the entropy of the 
flow line reaching the accretor with A^o = 1- The pressure 
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being uniform, Eq. (19) implies that the final Mach num- 
ber of the flow line indexed by w depends on its entropy 
S{'cu) according to: 



=4e«(^^i-s(-)) -3 



(69) 



The erroneous conclusion of Hunt (1971), that the axisym- 
metric flow with 7 = 5/3 never becomes supersonic close 
to the accretor is in contradiction with the 3-D numer- 
ical simulations (see Fig. 17 in Ruffert 1994b). It is even 
in contradiction with the Mach number J^o{9 = tt) ~ 2.1 
that one can deduce from the velocity, entropy and density 
displayed in Figs. 5, 6 and 8 in H unt (1971) for Moo = 2.4. 
In a similar way as Eq. (CIO), obtained for 7 < 5/3, 
Eq. (|9|) shows that the functions S, (3, M also cease to be 
independent at the boundary of the accretor for 7 = 5/3 
if — » 0. For a given distribution of entropy <S'(\E'), the 
only degree of freedom left by Eqs. (38), (67) and ( |69|) 
is the value of Si, i.e. the amount of mass accreted sub- 
sonically. The value of iSiopt which maximizes the total 
mass accretion rate satisfies an integral equation depend- 
ing only on S'(^'), which can be solved numerically. Thus 
if the shape of the shock is known, the distribution S{'^) 
is fixed by the Rankine-Hugoniot jump conditions along 
the shock (see Appendix^), and the total mass accre- 
tion rate M{S{'^), Si] is bounded by this maximum value 
M{5(«'),5iopt}. 

Since the subsonic accretion of material with a non- 
uniform entropy is possible for 7 = 5/3, we only have 
the inequality Ogo > dc- Nevertheless, the difference of en- 
tropies in the region of subsonic accretion is limited by 
the Eq. iM) to: 



Si-Siw = 0)<^log^ 



0.72 . 



(70) 



An angular sector of passing flow lines may exist, and is 



described in Sect. 4.2 



An estimate of the first order departures from spherical 
symmetry is obtained by writing the differential Eq. (27) 
to first order in the vicinity of r = (see Appendix ^ . 
This provides us with a lower bound on the departure 
from sphericity for the function /3, when the entropy is 
not uniform: 



< lim 



91og/3 
91ogr 



< 1 ff S'{^) / 



(71) 



As in Sect. 1.2, this "bending" of the flow lines enforced 
by the entropy gradient could be related to the forma- 
tion of long-lived radial shocks in the supersonic region of 
flows with a non-uniform entropy, as observed in numeri- 
cal simulations (Ruffert 1994b). 

5. Mass accretion rate for a spherically symmet- 
ric flow with non vanishing kinetic and thermal 
energies at infinity 



5.1. Analytical estimates of the mass accretion rate 

Only two situations are known to have analytical solu- 
tions: 

(i) Spherically symmetric accretion from a gas at rest 
at infinity (Bondi 1952). The mass accretion rate of the 
unique transsonic solution is: 



Mr 



AttG^AP 



37 



5-3t 
2(7-1) 



(72) 



(ii) Axisymmetric accretion from a gas with negligible 
temperature (Hoyle and Lyttleton 1939): 



Ah 



'i'nG^APr 



HL 



(73) 



Other normalizations were also proposed to account for 
the pressure effects with more realism, e.g. a factor 2 lower 
(Bondi & Hoyle 1944). 

In order to bridge the analytical gap between the case of 
spherical accretion described by Mb (Eq. ^2|) and the case 
of supersonic accretion from a gas of negligible tempera- 
ture described by Mhl (Eq. [t^), Bondi (1952) proposed a 
simple interpolation formula, meant to give the correct or- 
der of magnitude. Here we use the normalization adopted 
later by Shima et al. (1985) and Ruffert & Arnett (1994) 
on the basis of numerical simulations, which is a factor 2 
larger: 



M, 



BH 



(74) 



It has the advantage of being very simple and giving a 
roughly correct limiting behaviour for Voa — > or 
00. However, because it is independent of 7, it is bound 
to fail to reproduce quantitatively the analytical spherical 
accretion rate Mb, which varies by a factor ~ 3 between 
7 = 5/3 and 7 = 4/3. Moreover, numerical simulations 
have shown typical discrepancies of a factor 2-3 with this 
formula for various Mach numbers (see Fig. ^ and Fig. 
and also Hunt 1979). 

5.2. The mass accretion rate of a shocked spherical flow 

An intermediate situation, also analytically tractable, is 
the case of a spherical accretion with a constant kinetic 
energy at infinity: the gas is "thrown" towards the accre- 
tor, from a reservoir at infinity where the flow speed is 
Woo, the density p^o and the sound speed Cqo are uniform 
and non- vanishing, i.e. the kinetic energy and the thermal 
energies are constant at infinity. 

To this end, we imagine converging radial "pipes" of 
constant cross section (Fig. ^). They remain separate up 
to a minimum radius rp at which they merge. For a given 
distance rp, if a stationary solution exists, the mass accre- 
tion rate imposed at infinity would be: 



Al = Awr^pa 



(75) 
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_2^^-is^+AS)j^^ . (79) 

The mass accretion rate is clearly an increasing function 
of energy B^o, and a decreasing function of the total en- 
tropy (Soo + A5). 

We now consider the two possible cases of flows with and 
without a shock separately. 



5.2.1. Flows without a shock 



Fig. 6. The gas is "thrown" towards the accretor in a spher- 
ically symmetric way. "Pipes" of constant cross section are 
used to allow for both a constant thermal and kinetic energy 
at infinity. The accretion radius rA (full circle) can be smaller 
than the sonic radius ro (dotted circle). A spherical shock is 
represented by the thick dotted circle. 

Thus 2rp is also the radius of a cylinder containing the 
accreted material in an axisymmetric flow where the tran- 
sition between the parallel flow at inflnity and the spheri- 
cal flow near the accretor (and also the decrease of angu- 
lar momentum along each flow line) would be artiflcially 
driven through such pipes. This situation allows us to un- 
derstand how the kinetic energy influences the mass ac- 
cretion rate, for a gas with an adiabatic index 7. 
For spherical flows, the mass accretion rate of the station- 
ary solution is determined by the regularity condition at 
the sonic radius. We obtain exactly the same dimension- 
less quartic equation as in the classic spherical Bondi ac- 
cretion, determining, in the notations of Theuns & David 
(1992), r as a function of (see for example their Eq. 2.9), 
except that the velocities v, c and distance r are now nor- 
malized differently: 

1/2 



1 



7 ■ 



1 



-Mi 



GM 



-y— 1 r, 

7; 

2 00 



(76) 
(77) 



For spherically symmetric flows, the regularity condition 
at the sonic radius flxes the mass accretion rate Afgph for 
the unique transsonic solution: by integrating the mass 



flux ( p9[ ) over the sonic sphere (A^(ro) = 1), we obtain: 



sph 



"1' 


1+1 

2(7-1) 


[4(7 


-1)1 


2 




5- 


37 



2(7-1) 



AS r)2(7-l) 



(78) 



where AS is the entropy gained when passing through a 
spherical shock, if the flow at infinity is supersonic. 
Introducing the entropy 6*00 of the gas at infinity defined 
by Eq. (^), the spherical mass accretion rate can also be 
written: 



'1 


7+1 

2(7-1) 


[4(7 


-1)1 


5-37 

2(7-1) 
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37 J 





According to Eq. (78) when AS* = 0, we expect the mass 
accretion rate of isentropic flows to increase with the Mach 
number: 



M,ph(A5 = 0) = Mb 



1 



7-1 



Mi 



(80) 



This increase of the mass accretion rate with the Mach 
number for spherical subsonic flows, is not obvious a pri- 
ori. An increase of the speed at infinity acts as an addi- 
tional external pressure on the flow: according to Eq. (^3|), 
the sonic radius rp decreases, and Eq. (|29| ) shows that the 
momentum increases. However, the decrease of the area 
of the sonic surface, over which the momentum must be 
integrated to obtain the mass accretion rate, is more than 
compensated by the increase of momentum, for 7 < 5/3. 
The mass accretion rate is independent of the Mach num- 
ber for 7 — 5/3. 

5.2.2. Flows with a shock 

Increasing the velocity of the flow at infinity increases its 
total energy i?oo, but passing through a shock also in- 
creases its entropy. According to Eq. (79), these two ef- 
fects act in opposite directions on the mass accretion rate. 
Whether the mass accretion is eventually increased or de- 
creased depends on the position of the shock, which can 
be determined numerically as follows: for high Mach num- 
bers, we choose the distance of the pipes, so that the sound 
velocity prior to the shock ci equals the sound velocity at 
infinity ci = Cqo • This corresponds to the position rp > rgh 
such that the dilution due to gravitational acceleration is 
balanced by the concentration of flow lines. We make this 
particular choice by analogy with the case of axisymmet- 
ric accretion with a detached bow shock: we show in Ap- 
pendix^ that the density enhancement, along the axis of 
symmetry, before the shock, is smaller than a factor 2 if 
Tsh > 0.03rA. 

For each radius r, the Bernoulli equation (|^) determines 
the supersonic speed vi{r) and Mach number Mi{r) prior 
to the shock. Using the Rankine-Hugoniot jump condi- 
tions, we can compute the entropy jump AS{r) and the 
corresponding subsonic velocity V2{r) after the shock, if 
the shock occurs at that radius r: 



vi{r) 



2GM\ 



1/2 



r / 



Mi{r)cc 



(81) 



2 + {^-l)Mi , , 
V2ir) = —, — , -.N . .2 ^H'") 



-AS(r) 



(7 + l)-M? 

(7 + l)-M^ 
2 + (7-l)Al2 
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7 + 1 



2jMi -7 + 1 



(83) 



For a given entropy AS and Bernoulli constant Boo, we 
denote the velocity profile of the unique transsonic so- 
lution, by v[AS, Boo]ir). The unique shock position Tgh 
corresponds to the solution of the equation 



V2{rsh) = i'[AS'(rsh),Soo](r-sh) 



(84) 



Fig. 7. Position of the shock rah, in accretion radii, as a func- 
tion of yVfoo- The value of 7 is indicated on each curve. 



The position of the shock determined by this procedure 
is shown in Fig. ^, as a function of 7 and Adoo- Except 
for values of 7 close to 5/3, the spherical shock distance 
is larger than the accretion radius, and thus larger than 
was observed in numerical simulations of the axisymmet- 
ric BHL flow (typically Vgh ~ 0.2rA)- 
The spherical mass accretion rate is given by the Eq. (80) 
diminished by a factor exp(— AS'(rsh)) < 1 defined in 
Eq. (83). 

Increasing the kinetic energy of the supersonic flow al- 
ways decreases ultimately the mass accretion rate below 
the value AI-q obtained for a gas at rest. Between the two 
opposite tendencies (entropy and energy), the decrease 
due to entropy will dominate for large Mach numbers, 
since the effective Mach number A^i at the shock radius 
is larger than the Mach number at infinity. For Aioo ^ 1, 
Eqs. (79) and (83) imply: 



sph 



C(7) 



Mc 



Ml 



Pea ^ ri( \ 

-3- < C" 7 H~ 

via ^lo 



(85) 



Fig. 8. The mass accretion rate, as a function of Mao, in Bondi 
units Mb (upper curves) and also in units of the interpolation 
formula Mbh (lower curves). The value of 7 is indicated on each 
curve. The mass accretion rate follows the Hoyle-Lyttleton be- 
haviour when the kinetic energy dominates the thermal energy 
(slope A^""^). For a nearly isothermal flow, however, it exceeds 
the Bondi value if the thermal energy is larger than the kinetic 
energy. The dotted horizontal line is the mass accretion rate 
from a gas at rest at infinity. 



Moreover, Fig. |^ and Fig. || show that the shock radius is 
proportional to the accretion radius for high Mach num- 
bers, in other words the ratio Moo/ M.i is a non-vanishing 
constant for high Mach numbers. Apart from a multiplica- 
tive factor (7(7) depending on 7, the mass accretion rate in 
this spherical geometry scales like the mass accretion rate 
Mhl derived by Hoyle & Lyttleton (1939) for an axisym- 
metric flow of negligible pressure. It is remarkable that the 
mass accretion rate is asymptotically independent of the 
value of the sound speed at infinity. 
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5.3. What is the maximum spherical mass accretion rate ? 

We see from Eq. (79) that Mb is the maximum spherical 
mass accretion rate for 7 = 5/3. As we shall proceed to 
show, the situation is very different for 7 close to 1. The 
value of 7 defines the strength of a pressure change in 
response to a density change. Since the pressure gradient is 
the only force limiting the mass accretion rate, the higher 
the value of 7, the higher the pressure forces and thus the 
slower the accretion (see Eq. [72[ ). 

An isothermal gas is much more easily "swallowed" by 
the accretor: a significant increase of the mass accretion 
rate occurs for 1 < < 2/(7 — 1). In this range, the 

thermal energy dominates the kinetic energy, and Eq. ( p3[ ) 
implies that the sonic radius scales like the Bondi radius ^ 
GM/c^. More accurately, taking the Rankine-Hugoniot 
jump condition (83) in the limit 7 — > 1 and combining 
with Eq. (79) gives: 

M,ph(7 - 1) = MbMI exp ' + ■ (86) 

We can estimate the contribution of the exponential factor 
using the Bernoulh equation (81), for 7 1: 



2 rsh 
Consequently, the mass accretion rate scales like: 



(87) 



We conclude that for values of 7 approaching 1, the largest 
spherical mass accretion rate is: 



7 



1 



Mb 



(91) 



and is reached when the kinetic and thermal energies at 
infinity are comparable, i.e. : 



Mi 



1 ■ 



1 



(92) 



This scaling obtained in the isothermal limit for high Mach 
numbers is confirmed by the exact calculation of the mass 
accretion rate, for any value of 7, through the numerical 
derivation of the position of the shock, as shown by Fig. ||. 
The entropy-energy approach for a spherical flow has 
shown us that the mass accretion rate ultimately decreases 
like the Hoyle-Lyttleton rate when the kinetic energy 
dominates the thermal energy. 

Nevertheless, the mass accretion rate increases consider- 
ably above its value at M.00 = 0, when the kinetic en- 
ergy is comparable to the thermal energy, i.e. 1 < ~ 
2/(7 — 1). This effect is most pronounced for 7 ^ 1. Thus 
the value of the adiabatic index plays a very important 
role in the determination of the mass accretion rate in the 
case of spherical accretion with a shock. Similarly, we ex- 
pect the mass accretion rate of an axisymmetric BHL flow 
with constant parallel velocities at inflnity to depend on 
the adiabatic index. 



Msph(7 ^ 1) = 0{MbMI,) for 1 « Ml, < 



7 ■ 



1 



Using the geometrical constraint that rsh > ro, we can 
also write the mass accretion rate Mgh as the flux of mass 
through the shock surface, for 7^1. For high Mach num- 
bers, Eq. (81) ensures that vi ~ Uoo, and the density be- 
fore the shock pi ^ poo because the pipes are of constant 
cross section. Thus, 



,v _ . 2 AttG^M^p^ ( rsh 
Msh = 47rrsijPiWi = ^ 



Mc 



(89) 



we conclude from 
1/2 



Comparing this equation to Eq. ( |8q ) 

Msh = A^sph that the shock radius is a factor M__ 
larger than the sonic radius. Inserting this result back into 
Eq. (^, we further conclude that the mass accretion rate 
is: 



M,ph(7 ^ 1) = e 



_ 3/2 PooG^M^ 



-Mi 



for 1 < Mt. < 



7-1 



(90) 



which is a factor M^, larger than the spherical mass ac- 
cretion rate Mb from a gas at rest at inflnity, and a factor 
higher than the Hoyle-Lyttleton accretion rate Mbl- 



6. An interpolation formula for the mass accretion 
rate of an axisymmetric flow with 7 > 9/7 

The axisymmetric mass accretion rate Max can be for- 
mally calculated from the integral of the mass flux along 
any surface surrounding the accretor. We can decompose 
the mass accretion into three components, 



Afax = M,, 



Ms + My 



(93) 



where Mso is the mass flux through the sonic surface, Ms 
is the mass flux in the subsonic region (Types SF and 
FSF), and Mp is the mass flux in the supersonic region 
ahead of the shock (Type FSF only). 
The trajectories of the gas in the supersonic region ahead 
of the accretor are well approximated by hyperbolae, ob- 
tained by neglecting the temperature of the gas. M-p de- 
creases to zero when the accretor size r^, decreases to zero 
(Eddington 1926): 



Mp ^ Trr^irA + ?'*)Po 



(94) 



According to Sect. 4.4, Ms also decreases to zero with 



decreasing size of the accretor, for shocked flows with 7 < 
5/3. Consequently, the mass flux through the sonic surface 
gives a good estimate of the total mass accretion rate for 
small accretors: 



Maoc = Mso(l + e(r.)) for 1< 7 < 5/3 . 



(95) 
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Using Eq. (^9|) with = 1, we can write Mso as follows: 

1)" 



2(7- 
7 + 1 



5-37 



B' 



-AS 



(96) 



As in the case of spherical accretion, the mass accretion 
rate is an increasing function of the energy and a decreas- 
ing function of the entropy. However, the geometry of the 
sonic surface also influences the mass accretion rate. We 
know from Eq. (^2|) that the shape of the sonic surface 
rso{r,9) is directly related to the shape of the flow lines. 
In contrast to the energy and the entropy, the geometrical 
factor depends on the shape of the flow lines at infinity, 
and is more difficult to estimate since it does not corre- 
spond to a conserved quantity. One can try to obtain, for 
the axisymmetric case, an interpolation formula for the 
mass accretion rate Mj, based on the spherical entropy- 
energy approach. 

With some trigonometry we write: 



d$ 



cos (3 — sin (3 



d\ogrso\ sin 9 



-d0 



(97) 



47rr^ 

' so 

The direction 6 — being in the direction of the flow, the 
axisymmetric sonic surface is described by rso{d), with 
< 6' < For Type F flows, 6*50 must be replaced by tt 
in what follows: 



An 





7+1 


'I 


2(7-1) 


2 





1 



2 '00 
5-37 



h{f) ^ ^ • ' . (99) 

First note by differentiating h{f) with respect to r that: 
h{f) > h{fo) = 1 , (100) 

The integral in Eq. (98) can be seen as an average of /i(fso) 
times the entropy weighted by the geometrical function 
d$/47rrgQ > 0, thus defining an effective entropy corre- 
sponding to an effective Mach number A^o = ■Mi{9c)- 

d$ 



1 



d$ 

47rr2^ ' 

7+1 

2(7-1) 



(98) 



47rrL 



d$ 



(101) 



Using Eq. (83) and the definition (||) oi Mb, Eq. (98) be- 



comes: 



2 + ij-l)Ml 



7 + 1 
27>12 _ ^ + 1 



d$ 



47rr2 



.(102) 



One can check that for isentropic flows (AS* = 0) and 
small deviations from sphericity {f{6) ^ f3 ^ {fgo — fo) <C 
1), the mass accretion rate is accurately described by the 
spherical model of Sect. ||: 



M,o = A4ph(A^oo) [1 + {Max(f,o - ro)'}] 



(103) 



In order to approximate the mass accretion rate Mgo, one 
must estimate the value of the effective Mach number 
A^o, the effective sonic radius fso(^'e) and the integral of 
{d^ / 4:7rr^^) . If the shock is detached, we can make simple 
guesses if the spherical sonic radius tq is smaller than the 
accretion radius ta. We then use the spherical approxima- 
tion for $ and rgo, and use an ad- hoc expression for A4o 
only, to account for the asymmetry. We see from Fig. ^ 
that ro < TA for 5/3 > 7 > 9/7. The range 9/7 > 7 > 1, 
and especially the vicinity of 7 = 1, requires a more care- 
ful treatment because ro can be orders of magnitude larger 
than ta. Consequently, we restrict the present approach 
to the range 5/3 > 7 > 9/7, and make the simple approx- 
imation fsoi^e) ^ ^Q. 

Since we would like to recover the Hoyle-Lyttleton accre- 
tion rate Mhl for high values of the Mach number, we see 
from Eq. ( |85|) that this requires a constant value for the ra- 
tio A^c/A^oo when A^oo — 00. As a first approximation, 
we propose the ad-hoc prescription that Mso AA/hl 
when Aloo — 00, and retain the same ratio Ade/Moa for 
intermediate values of A^oo- Neglecting the variations of 
the sound velocity before the shock, we can write the effec- 
tive Mach number in terms of the angle tpc of the "effective 
flow line" of entropy S{A4c) and the normal to the shock 
surface, and the distance between the accretor and the 
intersection of this flow line with the shock surface: 



Ma 
Aloe 



COS (fc 1 



(104) 



A geometrical interpretation is the following. Let the ac- 
cretion radius act as a scaling factor on the topology 
of the accreted flow, and keep the angle (pc and the ratio 
fK/'''c independent of the Mach number. Eq. (104) then 
implies that the ratio A4c/A^oo is constant. This does not 
contradict the fact that the opening angle of the Mach 
cone, well after the accretor, depends on the Mach num- 
ber, since our approximation deals with the region of the 
bow shock which is mainly ahead of the accretor. 
We obtain, for Mc > 1: 



10 

rso{9e 

1 



Mac 2iX^ 



d$ 
47rr2^ 

)^fo 

2" 

7_ 



(7+1)- 



(7-1)^(5-37)^ 



(105) 
(106) 

(107) 



With A^c defined by this relation if Aic > 1, and 
Mc = 1 otherwise, our interpolation formula follows from 
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Eq. (102): 



-Mi 



_2 + {j-l)Ml 



1 



1 



27^1-7 + 1 



(108) 



Fig. 10. Same as Fig. |, but for 7 = 4/3. 



Fig. 9. Comparison between our interpolation formula and the 
numerical simulations. Mass accretion rates are in A/bh units. 
The solid line corresponds to A = 1, the dotted line to A = 0.5. 
For 7 = 5/3, the numerical simulations provide only upper 
bounds to the mass accretion rate of a point mass because the 
sonic surface is always attached to the accretor. 



We extend our rough approximation to the case 7 = 
5/3, neglecting the mass accreted in the subsonic region. 
We know from Sect. 4.6 that this corresponds to a small 



(but not zero) fraction of the total flux of mass. 
The curves displayed in Fig. ^ and Fig. 10 show the mass 
accretion rates obtained with our interpolation formula 
Mi, the Bondi-Hoyle interpolation formula Mbh and 
the results of numerical simulations (Ruffert 1994b, 1995; 
Ruffert & Arnett 1994). For 7 = 5/3 our interpolation 
Eq. (108), in Mbh units, rises monotonously with Mach 
number (Fig. The numerical results fit only marginally 
well. Note that the numerical simulations with 7 = 5/3 
give only an upper bound for the mass accretion rate, 
since the accretor is always larger than the spherical sonic 
radius. For 7 = 4/3 (Fig. a prominent maximum is 
apparent both on the interpolation formula and on the nu- 
merical results. Comparing Fig. |l^ with the corresponding 
curve in the lower part of Fig. |^ gives an indication on the 
role played by non-sphericity in reducing the mass accre- 
tion rate. 

Despite the considerable simplification of ignoring the ge- 



ometrical term $(0) and keeping the ratio A^c/A^oo inde- 
pendent of A4oo, these curves are in reasonable agreement 
with the numerical simulations for 7 = 4/3 and 7 = 5/3. 
Thus we expect this analytic formula to be relevant for 
the whole range 5/3 < 7 < 9/7. 

Note that we do not solve the question of the dependence 
of the mass accretion rate with respect to the adiabatic 
index for high Mach numbers. For this purpose, and fol- 
lowing Bondi & Hoyle (1944), we introduced the constant 
A, assumed to be between 0.5 and 1. Our interpolation 
formula is nothing more than a natural way to link the 
Bondi mass accretion rate at A^oo =0 with the Hoyle- 
Lyttleton mass accretion rate at M.00 ~* 00. 
One should also keep in mind that the mass accretion rate 
obtained in the numerical simulations is for most models 
a time average over the fluctuations or quasi periodic os- 
cillations due to the instability of the flow, and therefore 
differs from the exact mass accretion rate of the solution 
of the stationary equations. This instability is stronger for 
7 close to 5/3. 

7. Conclusions 

By using a geometrical interpretation of the shape of the 
flux tube at the sonic radius, we have shown the partic- 
ular signiflcance of the radius tq (Eq. ^3|) , which depends 
only on the adiabatic index and the energy of the flow. 
It is not only the sonic radius of all spherical flows, but 
is also related through the property V (Sect. |3.3D to the 
sonic surface of axisymmetric flows. 

By reducing the stationary flow equations to a single par- 
tial differential equation, we have performed a local anal- 
ysis near the accretor and have obtained some general 
results which aim at clarifying the diversity of possible 
configurations. We have extended the classification intro- 
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duced by Bondi (1952) to axisymmetric stationary flows. 
In addition to the subsonic and supersonic types of accre- 
tion, angular sectors without any accretion are also pos- 
sible, as well as singular directions of infinitely high ac- 
cretion flux. This is illustrated by the self-similar solution 
found by Bisnovatyi-Kogan et al. (1979), which is a very 
particular singular case {9c — 0). 

The angular sectors without any accretion is subsonic for 
flows with a detached shock, so that 6so < 9c in the nota- 
tion introduced in Sect. 4.1. 



Among some of the properties that we established for flows 
with 7 < 5/3 are the following: 

(i) the sonic surface is likely to be attached to the 
accretor for nearly isothermal flows with high Mach num- 
bers, 

(ii) isentropic accretion proceeds mostly from the 
downstream hemisphere of the accretor, 

(iii) subsonic accretion must be isentropic. 

We have also found a series of properties for flows with 
7 = 5/3: 

(i) accretion is always regular, 

(ii) the sonic surface is always attached to the accretor, 

(iii) the pressure distribution is spherically symmetric 
to first order, 

(iv) the Mach number and the entropy at a point like 
accretor are algebraically related to the local mass flux, 

(v) shocked matter is accreted mostly from the down- 
stream hemisphere, where the density, velocity and Mach 
number are maximum, and the temperature and entropy 
are minimum, 

(vi) isentropic accretion is spherically symmetric to 
first order for the velocity, temperature, density, mass flux, 
and Mach number. 



noted that the mass accretion rate scales like the Hoyle- 
Lyttleton formula for high Mach numbers. Moreover, we 
have outlined the particularity of nearly isothermal flows 
for which the mass accretion rate can be a factor 2/(7—1) 
higher than the Bondi mass accretion rate, when the ki- 
netic energy at infinity is comparable to the thermal en- 
ergy. This demonstrates the fundamental difference of the 
BHL flows that have 7 close to 1 and those with 7 close 
to 5/3. 

We have shown that the total mass accretion rate equals 
the mass accretion rate through the sonic surface, for 
shocked regular flows with 7 < 5/3. We have derived an 
interpolation formula for the mass accretion rate valid in 
the range 9/7 < 7 < 5/3. It links spherical accretion to 
the case of highly supersonic accretion, with reasonable 
agreement with numerical simulations. 

Analytical estimates of the following quantities, how- 
ever, are still missing: 

(i) the distance of the shock as a function of 7 and 

Moo, 

(ii) the value of 6so and 6c as a function of 7 and Aioc, 

(iii) the mass accretion rate in the limit of infinite 
Mach numbers, as a function of 7. 

Nevertheless, our classification will help perform, in 
a forthcoming paper, a local stability analysis of all the 
configurations considered. 
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3-D numerical simulations have shown that the 
strongest instabilities occur with the adiabatic index 7 = 
5/3. The fact that subsonic regions invariably reach the 
accretor for 7 = 5/3 might be considered as a warning 
about the possible role, in the instability mechanism, of 
the boundary conditions at the surface of the accretor. For 
7 < 5/3 and sufhciently small accretors, it would be inter- 
esting to check for the presence of the instability in cases 
where the accretor is surrounded by a su personic region 
(Type F in the nomenclature of Sect. 3.1). 



We have also determined the departure from spheric- 
ity, and constrained the radial variation of the angle (3 for 
r ^ 0. In particular, we have stressed the influence of the 
entropy gradient on the bending of the flow lines for su- 
personic flows with 7 < 5/3 and 7 ~ 5/3. 

We have stressed the role of the geometrical shape of 
the shock in the determination of the entropy of the ac- 
creted matter, and contrasted it with the simpler case 
of shocked spherical flows. In this latter case, we have 
shown the relevance of the entropy-energy approach, and 



A. Third order partial differential system for an 
ELxisymmetric stationary flow 

The combined vorticity equation, continuity equation and en- 
tropy equation are written as functions of A^, /3, 5 as follows: 



d\ogr 



aO 



2 + (7 - l)M'^ de 
sin^ /3 sin /3 cos jS 
tan e 1+f 

dlogM 



+ 



l^M^ dS 



+ 



'jM^ de 
2-7 1 + 2f 



7-1 2 

2 M'^ sin [3 cos 13 dlogM^ 



(Al) 



2 + (7-l)7W2 aiogr 2 + (7-l)7W2 
5-37-2(2-7) sin^/3 



dp tan f3 dS 
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-(2cos^/3 + 1) 
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dS 
d log r 



2(l + f)(7-l) 
sin (3 cos j3 
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= — tan/3 — 



(A2) 
(A3) 
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B. Proof of the property of the axisymmetric 
sonic surface of Type F flows 

Using the axisymmetry condition and some simple trigonome- 
try, the variation of the cross-section $ along the flow line can 
be written: 



d log $ dl 
dl d log r 



tan B 
tan f 



dp 



(Bl) 



cos (5 dn ' 

where j3{r, 9) is defined by Eq. (|ll|), and its derivative is taken 
in the direction n perpendicular to the flow line (oriented to- 
wards increasing 6). For P — we recover the spherical case 
of conical flow tubes, with a cross section proportional to the 
square of the distance. 

We can use Eq. (Bl) along the detached sonic surface, together 
with Eq. (|3|), for < 6* < tt: 



m = ^ 



^0 - fW 

1 f Tso dp 



with 



+ 



tan P 



2 y cos P dn tan 9 



(B2) 
(B3) 



Let us now prove that the sonic radius reaches the value of the 
spherical radius at least once, i.e. that f{9) has to vanish. 
Axisymmetry imposes P{r,9 = 0) = P{r,9 = -k) = Q, and by 
Eq. (B3): 



/(O) 



(B4) 
(B5) 



Suppose, for example, that the sonic surface is inside the sphere 
of radius ro and detached from the accretor, i.e. f{9) > e > 
for 9 G [0,7r]. Equations (B4) and (B5) imply that P{r,9) is 
strictly positive near Tsoifl — 0), and strictly negative near 
rso{9 = tt). A line must exist in the {r,9) plane separating 
rao(O) from rsoiir), along which P vanishes, when decreasing 
from a region of positive values to the region of negative values. 
We denote by rso{9i) the intersection of this line with the sonic 
surface: 



dp 



cos P dn 



< 



(B6) 



which contradicts the hypothesis. The same contradiction 
arises if we suppose that f{9) < —e < 0. Consequently / must 
vanish for some value of 9. 

C. Linearized partial differential equation for 7 < 
5/3 

By using Eq. (23) , we can relate the second radial derivative of 
^ to its first radial derivative, at r = 0, through the following 
equation: 



lim 









— lim 




or 



aiog/3 _ ^ ^ _ (^\ 



9 log 



dfd9 \d9 ) 



(CI) 



The local mass flux is finite for regular flows, and therefore the 
limit when r — > of (fd^'^ /dfdff) is zero in the angular sectors 
where the flow is regular. 



Let us define the new radial variable z 
tion (129) implies to first order: 
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( 1 d-fo \ 
{sm9 d9 J 
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(i) If Eq. ( jSq ) is a strict inequality, Eq. (^) implies that the 
pressure is spherically symmetric to first order. Using Eqs. (38) 
and (52), the stream function 'to(S) at the surface of the ac- 
cretor is directly related to the entropy So{9) by: 
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The leading terms of the differential Eq. (27) satisfy the fol- 
lowing equation: 

d^ d 
(1 + + (1 + ^2)b— + (1 + e3)c + 

+ + + ^^)^} *^ = + ' 

where the six functions ei{z,9) converge to zero when 2^0. 
The functions a, b, c, d depend on 9 only, and are defined as 
follows: 
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(C6) 
(C7) 
(C8) 



Assuming lim2^o(91og ^'i/91og 2) > 2 and d 7^ contradicts 
Eq. (C4). Thus we conclude that: 



91og*i 
d log 2 



< 2 for ^'(^o) / 



(C9) 



By derivating Eq. (27) with respect to r and considering 
the leading terms, we find that zd'^\/dz satisfies the same 
Eq. (C4) as Vpi does. Using Eq. (23), we obtain Eq. (jH^ for 
regular sectors with a non-uniform entropy. 



(ii) If Eq. ( |56[ ) is an equality, we use Eq. (CI) in order to 
write the first order terms of Eq. (27) as follows: 



,. 9* 

z^Q dz 



4e^o sin 9 



d 



37(5-37)(2-7) d9 



Vsmi 
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In both cases (i) and (ii), So,'^o and the radial derivative 
[d"^ /dz)o cease to be a set of three free functions of 9 on the 
surface of the accretor, when the accretor size decreases to zero. 
Let us recall that for spherical flows, one parameter (the mass 
flux) is required as boundary condition on the surface r = r* 
of the accretor, with no particular change if 0. In this 

respect, the differential system (27) is more singular at r = 
for axisymmetric flows than for spherical flows. 
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D. Linearized partial differential equation for 7 = 
5/3 

Because the stream function \[' is finite everywhere, e* and 
f^a^^/^f^ tend to zero with decreasing r. Eq. (^ ) implies 
that when 7 = 5/3, the Mach number tends to a finite limit 
everywhere where the local mass flux (i.e. d'^o/dO according to 
Eq. (38)) is not zero. Taking the limit of Eq. (27) when f — > 0, 
we find that its last term must vanish at f = 0: 



1 9^0 
ta.ne~de' 



96(2 V de 



Ml 



0,(D1) 



where the subscript denotes the value of the function at r = 0. 
This equation is equivalent to the condition that the pressure 
is spherically symmetric to first order. 

In order to find the behaviour of the quantities ^',A^^ near 
f = 0, we consider only the leading terms of Eqs. (^6|) and 
(27). We first define Mi and ^Pi as follows; 



M^{r 



= Ml{e){l + Mi{r,e)) 
*(r,6l) = *o(6') + *i(r,6l) . 



(D2) 
(D3) 



By linearizing Eq. ( pq ) near Ml 7^ 1, we obtain to first order 
(3 + Xg)* 9*0 a*i 



2Ml sin^ 



de 



(D4) 



Introducing this into Eq. (27) , we obtain the following leading 
terms; 

(1 + ^^y-Q02 + (1 + ^8)/^ + (1 + ^9)5 + 

(1 + eio)r'|^ + 1(1 + eiO^Ji 1 *i = (1 + ei2)ftf , (D5) 
ar^ I or \ 

where the six functions ej{f, 9) converge to zero when f ^ 0. 
The functions e, f,g,h depend on 9 only, and are defined as 
follows; 



1 — Mq tanf 



89 5Ml 



{3 + Mlf /9*o 
3MI{1 - Miy sin'' 9 \ 89 



^ (8<foY 3g"(*o) 
^ - \ 89 J 5Ml ' 

4 S + Ml / 9v&o \" g'(*o) 
~ 5{l~Miy\89) Ml • 

Let us suppose that I is a positive number such that 



lim 



8 log r 



(D6) 

(D7) 
(D8) 
(D9) 

(DIO) 



If the entropy is not uniform {h 7^ 0), the first two terms of 
Eq. (D5) are negligible compared to the rhs term. Dividing the 
remaining equation by f, we obtain; 



& 



8 



(1 + ^7)675^ + (1 + ^s)fj^ + (1 + eg) 



According to Eq. (DlC), lim,.^o*i/r = 0. Thus Eq. (Dll) 



contradicts our hypothesis (DIO), and we conclude that if 



< 



lim^|<l 

r^o dlogr 



(D12) 



By derivating Eq. (27) with respect to f, and keeping the lead- 
ing terms, fd'^i/8f satisfies exactly the same equation (D5) 
as ^'i does. From the same reasoning, we conclude that 



< lim 



9 log 8'i!i 
->(> 8 log r 8 log f 



< 1 



(D13) 



This lower bound on the scaling of f8'i!\/8f is also a lower 
bound on the scaling of f3: we obtain Eq. ( [7l| ) by using 
Eqs. (23), (D13) and the regularity of the fiow for 7 = 5/3. 



E. Approximation of the orbits within the super- 
sonic region by hyperbolae 

We check the validity of neglecting the pressure forces in the 
supersonic fiow ahead of the shock, and assuming that the 
sound velocity is uniformly constant to its value at infinity Coo • 
The velocity of the supersonic incoming gas is then defined 
by the unique hyperbolic trajectory passing trough this point, 
with velocity at infinity Vao ■ 

For each point (r, 9) of the shock surface, we can compute 
numerically the dimensionless impact parameter a (normalized 
to the accretion radius ta) of this hyperbola; 



2a^ 



l + (l + 4a2)i/2cos(6l- 



The half angle Baa of the hyperbola is defined by 



cos 9a 



(1 + 4a2)i/2 



(El) 



(E2) 



Let us calculate the variation of density along the axis of sym- 
metry, due to both the compression by gravitational conver- 
gence of the hyperbolic trajectories and the dilation due to ac- 
celeration. This can be done using the continuity equation (^ 
in cylindrical coordinates, along the axis, for 9 —> tt: 



8pf^v. 



(E3) 



The azimuthal derivative of the azimuthal component of the 
velocity is directly related to the impact parameter a through 
the conservation of angular momentum, while the (negative) 
radial component of the velocity, along the axis, is given by the 
conservation of energy: 



rve = —avc 

Vr{9 — Tv) — —Voo 



1 + 



1\ 2 



(E4) 
(E5) 



In order to determine the azimuthal dependence of the impact 
parameter, let us consider the trajectories defined by Eq. (El) 
for (9 — tt) ^ 0. Since a(7r) = 0, we write; 



{l+ei2)h .(Dll^(^) ^ _ n)^(n) + O ({9 ^ n)^) 



(E6) 



20 



An analytic study of Bondi-Hoyle-Lyttleton accretion. I 



From first order expansions of Eq. @ and Eq. @, we find, 
for (e - tt) ^ 0: 



(day 



^ da 



(E7) 



Choosing the negative solution of this equation (a is positive 
and decreases to zero when 6 increases to tt), leads to: 



dve , Va 



1+1 



f \ 2 



Using Eqs. (g8|) and (E5) in Eq. (|E^) we obtain: 
aiogp , 



df 



1 f + 2r 



(E8) 



(E9) 



This can be integrated to obtain the radial dependence of the 
density along the axis of symmetry: 



(ElO) 



The compression effect due to the convergence of the flow lines 
is consequently stronger than the dilatation due to the accel- 
eration of the flow. 

If the shock is not too close to the accretor, gravitational focus- 
ing is negligible and the density and temperature can be ap- 
proximated as consta nt in the upstream supersonic flow. One 
can check from Eq. (ElO) that this density enhancement is 
smaller than a factor of two for f > 0.03. Note that when the 
shock is detached, its distance is generally of order£sh(7r) ~ 0.2 
in numerical simulations. We deduce from Eq. (ElO) that the 
sound velocity scales like f''^"^-'''^ when f ^ 0, and is there- 
fore always negligible compared to the gravitational potential 
scaling like f'^. 

F. Entropy gradient produced by a detached shock 
in the vicinity of the axis of symmetry 

Let (a;ah,2/sh) be local coordinates in a frame chosen in a way 
such that Xsh is perpendicular to the shock surface, and ysh is 
parallel to it. Let us call A4i the Mach number relative to the 
component of the velocity, ahead of and perpendicular to the 
shock surface. In our local system of coordinates. 



Ml 



Cl 



(Fl) 



We shall use the classical convention of using the index "1" for 
quantities ahead of the shock, and the index "2" for quantities 
just after the shock. 

We write the entropy gradient after the shock as: 



V2 882 _ V2 882 dMi 

Vx2 dy Vx2 dMi dy 



(F2) 



We can use the classical jump conditions (83) to express the 
entropy gradient as a function of TVfi only: 



v{Mi 



VS = where 
Vx2 J —1 oy 

27(7 -1)(A1?- 1)2 



(2 + (7-l)>lf)(27>l?-(7-l)) 



< 1 



(F3) 
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We obtain Vxi by the vector product of the velocity and 
the vector tangent to the shock surface. Let the supersonic 
flow before the shock be described in cylindrical coordinates 
by the flow velocity Vr{r, 6), veif, 6). For the sake of clarity, we 
denote the velocity along the shock surface rah(S) by v{6) = 
v{rsb{d), 6) and a dot stands for the derivative with respect to 
the angle 6: 



vxi{e) 



(F5) 



We now prove that a regular axisymmetric shock surface can- 
not produce a local region of uniform entropy behind it, and 
more precisely, that the gradient of entropy created by the 
shock scales like the square of the azimuthal angle {9 — n) in 
the vicinity of the axis, for any shape of the shock which is 
regular and symmetric. 

When not specified, the quantities are taken at 6^ = tt in what 
follows. We also define Vgh = Vr{n) and r = fsh(7r). The sym- 
metry condition imposes: 



rsh(7r) 







(F6) 



Using the symmetry of the problem, we make the following 
expansions for small values of {0 — n): 



Vr{9) = Vsh + 



{9- 



df 



+ 



9V 



+0{{9-^f) , 

Me) = {e-^)^ + o{{9-nf) , 
f.^{9) = + ^^^f.h + o ((e - ^)*) . 

Using the expansions (F7), (F8), and (F9), we obtain 
Vxi{9) 



(F7) 
(F8) 

(F9) 



1 - 



fsh \ ^ Tsh / 2 dvg fsh dVr 



+0{{9~n)') 



1 d'^Vr 

'^~d9^ 
(FIO) 



In order to evaluate this expression, we use the approximation 
of the supersonic flow Vr{f,9),vg{fj9) by hyperbolic trajecto- 
ries. Derivating the Bernoulli Eq. (X^ with respect to f and 9, 
we obtain for 9 = tt: 



89^ ^ ~«eh \89 ) 

dVr ^ via 

8f 2ushr2^ 



(FU) 
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Replacing these expressions, together with Eqs. (E5) and (E8) 
into Eq (FIO), the effective Mach number is: 



Mi{9) 



{9~.f 



3 + 4f,h + 4fJ^(l + r.h)f;^^2 



16fsh(l + rsh) 



+ 0{{9-nf) 



(F13) 
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where the positive contribution depends on the second 
derivative fsh.i'^) of the shock surface on the symmetry axis, 
and is defined as: 

^ _ rsh _|_ J_ ^dve_ _ rah dvr \ 

We conclude that for any shape of the detached shock surface, 
the effective Mach number near the symmetry axis decreases 
at least as the square of the azimuthal angle {6 — tt)^. 
This result is particularly useful in order to prove that the 
mciss accretion rate in the subsonic region of Type SF flows 
decreases to zero with the size of the accretor. 



(F14) 
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